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Direct numerical simulations (DNS) are used to examine the pressure fluctuations gen- 
erated by a Mach 6 turbulent boundary layer with nominal freestream Mach number of 6 
and Reynolds number of Re T ~ 464. The emphasis is on comparing the primarily vortical 
pressure signal at the wall with the acoustic freestream signal under higher Mach number 
conditions. Moreover, the Mach-number dependence of pressure signals is demonstrated 
by comparing the current results with those of a supersonic boundary layer at Mach 2.5 and 
Re T « 510. It is found that the freestream pressure intensity exhibits a strong Mach number 
dependence, irrespective of whether it is normalized by the mean wall shear stress or by 
the mean pressure, with the normalized fluctuation amplitude being significantly larger for 
the Mach 6 case. Spectral analysis shows that both the wall and freestream pressure fluc- 
tuations of the Mach 6 boundary layer have enhanced energy content at high frequencies, 
with the peak of the premultiplied frequency spectrum of freestream pressure fluctuations 
being at a frequency of luS/Uoo 3.1, which is more than twice the corresponding frequency 
in the Mach 2.5 case. The space-time correlations indicate that the pressure-carrying ed- 
dies for the higher Mach number case are of smaller size, less elongated in the spanwise 
direction, and convect with higher convection speeds relative to the Mach 2.5 case. The 
demonstrated Mach-number dependence of the pressure field, including radiation intensity, 
directionality, and convection speed, is consistent with the trend exhibited in experimental 
data and can be qualitatively explained by the notion of ‘eddy Mach wave’ radiation. 


Nomenclature 

C p heat capacity at constant pressure, J/(K-kg) 

C pp Space-time correlation coefficient of the pressure field, dimensionless 
C v heat capacity at constant volume, J/(K-kg) 

H shape factor, H = 6* /9, dimensionless 
M Mach number, dimensionless 

M r relative Mach number, M r = ({Too — U s )/aoo, dimensionless 
Pr Prandtl number, Pr = 0.71, dimensionless 
R ideal gas constant, R = 287, J/(K-kg) 

Reg Reynolds number based on momentum thickness and freestream viscosity, Reg = p ” (/ °° , dimensionless 
Res 2 Reynolds number based on momentum thickness and wall viscosity, Res 2 = , dimensionless 

Re T Reynolds number based on shear velocity and wall viscosity, Re T = Pw ™' rS , dimensionless 
T temperature, K 

T r recovery temperature, T r = T oa ( 1 + 0.9 * K 

Ub bulk convection speed of pressure-carrying eddies, m/s 
U s convection speed of ‘effective’ radiation sources, m/s 
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Uoo freestream velocity, m/s 

a speed of sound, m/s 

p pressure, Pa 

q dynamic pressure, Pa 

s entropy, J/K 

u streamwise velocity, m/s 

u T friction velocity, m/s 

v spanwise velocity, m/s 

w wall- normal velocity, m/s 

x streamwise direction of the right-hand Cartesian coordinate 
y spanwise direction of the right-hand Cartesian coordinate 
z wall-normal direction of the right-hand Cartesian coordinate 
z T viscous length, z T = v w /u T , m 
7 specific heat ratio, 7 = C p /C v , dimensionless 

6 boundary layer thickness, m 

S* displacement thickness, m 

k thermal conductivity, k = pC p /Pr , W/(m-K) 

9 momentum thickness, m 

p dynamic viscosity, p = 1.458 x 10 ~ 6 777 x^ 4 , kg/(m-s) 
v kinematic viscosity, v = p/ p, m 2 -s 
p density, kg/m 3 

ui frequency, rad/s 

vorticity, s _1 
Subscripts 

i inflow station for the domain of direct numerical simulations 

rms root mean square 

w wall variables 

00 freestream variables 

t stagnation quantities 

Superscripts 

+ inner wall units 

(•) averaged variables 

(•)' perturbation from averaged variable 


I. Introduction 

It is well known that the free-stream disturbances in a conventional high-speed facility are typically 
dominated by acoustic radiation from the turbulent tunnel wall boundary layers . 1-13 The elevated freestream 
disturbance levels in conventional (i.e., noisy) high-speed wind tunnels usually result in an earlier onset of 
transition relative to that in a flight environment or in a quiet tunnel. Yet, the conventional facilities continue 
to be used for transition sensitive measurements because of the size and Reynolds number limitations of 
existing quiet facilities and the prohibitive cost of flight tests. To enable a better use of transition data from 
the conventional facilities, it is important to understand the acoustic fluctuation field that dominates the 
freestream disturbance environment in those facilities. Measurements at Purdue University indicate that the 
effect of tunnel disturbances on transition cannot be quantified in terms of a single metric corresponding 
to the root-mean-square amplitude of freestream acoustic disturbances . 14 With increased knowledge of the 
receptivity mechanisms of high-speed boundary layers , 15, 16 it becomes particularly important to characterize 
the details of the tunnel acoustics originating from the tunnel-wall turbulent boundary layers. Even during 
flight, the laminar boundary layer flow over the aerodynamic surface of interest may be exposed to acoustic 
radiation from the turbulent boundary layer over other neighboring surfaces, not unlike a test model mounted 
within a conventional wind tunnel. 

With the exception of early measurements by Laufer 1 and a few others , 5 all previous studies of pressure 
fluctuations due to high-speed turbulent boundary layers have focused on the surface pressure fluctuations 
underneath the boundary layer, and the existing data exhibit a considerable degree of scatter . 17-19 As 
pointed out by several authors , 19,20 there are few (if any) reliable measurements of the variance of the wall 
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pressure fluctuations and its frequency spectra at high speeds, due to either a poor spatial resolution of the 
pressure transducers and/or the limited frequency response of these sensors. Laufer’s measurements of the 
acoustic fluctuations in the freestream region 1 are also subject to analogous sources of experimental error. 
Moreover, as noted by Laufer, the interpretation of disturbance measurements in a wind tunnel is further 
complicated by the fact that the measurements reflect the combined outcome of acoustic radiation from all 
sides of the tunnel wall. As a result, highly accurate measurements of the absolute amplitudes of the radiated 
acoustic energy were not pursued during his experiments and only the statistical quantities that were least 
likely to be influenced by the presence of multiple tunnel walls were investigated. 

Direct numerical simulations (DNS) can overcome many of the aforementioned difficulties with exper- 
imental measurements of freestream disturbances as well as provide access to quantities that cannot be 
measured easily. The DNS can also isolate freestream disturbances due to acoustic radiation, thereby avoid- 
ing any contamination due to secondary sources such as vortical and entropy fluctuations in the incoming 
stream. The successful application of DNS for capturing the freestream acoustic pressure fluctuations has 
been demonstrated by the current authors in the context of a simulation of a Mach 2.5 supersonic turbulent 
boundary layer, in which many single and multi-point statistics of the pressure field, including the pressure 
intensities, frequency spectra, space-time correlations, and convection speeds, were described. 21 The current 
study extends that previous study into the hypersonic Mach number regime (where even the properties of 
surface pressure fluctuations have not been quantified as yet) to highlight the variation of acoustic charac- 
teristics with Mach number as well as to provide additional validation of simulations against experimentally 
measured trends. 

The flow conditions selected for the numerical simulation of a hypersonic turbulent boundary layer and the 
numerical method used for this simulation are outlined in Section II. Section III is focused on an analysis of 
the statistical characteristics of freestream pressure fluctuations, including power spectral densities, pressure 
fluctuation intensities, integral length scales, and convection speeds. Conclusions thus far are outlined in 
Section IV. 


II. Flow conditions and numerical methodology 

Table 1 outlines the flow conditions for the Mach 6 simulation. These conditions are similar to the 
operation conditions of Purdue Mach 6 Quiet Tunnel. Flow conditions for the previous DNS simulation of 
Mach 2.5 boundary layer are also included to facilitate the comparison between different Mach number cases. 
Observe that the values of Re T and Res 2 are similar for both DNS cases, allowing a meaningful comparison 
between the two simulations. T w /T r is approximately 0.8 for the Mach 6 simulation and 1 for the Mach 2.5 
simulation. 

Figure 1 shows the general computational set-up for the DNS in the present work, which parallels the 
setup in Duan et al. 21 for the Mach 2.5 simulation, wherein the effects of domain size and grid resolution 
were also assessed. The choice of grid parameters for the present study is based on lessons learned from 
Duan et al., 21 allowing a single simulation to be used to assess the effects of Mach number. The domain 
size for the present computation is (58.7<5i, 15.7(5 i, 39.7(5Q in the stream-wise, span-wise, and wall-normal 
directions, respectively. The streamwise domain length, L x = 58.7(5^, is selected to be larger than the 
eddy decorrelation distance to guarantee minimal spurious correlation being introduced due to the inflow 
turbulence generation. The number of grid points is (1600, 800, 500) in the streamwise, spanwise, and wall- 
normal directions, respectively. Uniform grid spacings are used in the streamwise and spanwise directions 
with grid spacings Ax + = 10.2 and Ay + = 5.5, respectively. The grids in the wall-normal direction are 
clustered in the boundary layer with A z + = 0.55 at the wall, and kept uniform with A z + = 5.7 in the 
freestream until up to approximately 6Si. Such grids are designed to adequately resolve both the boundary 
layer and the near field of acoustic fluctuations radiated by the boundary layer. Unless otherwise stated, 
the grid resolutions specified in this section are normalized by the viscous length scale z T at the inflow 
plane. The computational grid resolution inside the boundary layer is comparable to those reported in the 
literature in the context of previous simulations of turbulent wall-bounded flows using comparable numerical 
algorithms. 22-25 

The details of the DNS methodology, including numerical methods, initial and boundary conditions, 
has been documented in our previous paper. 21 Therefore, only a cursory description is given here. The 
full three-dimensional compressible Navier-Stokes equations in conservation form are solved in generalized 
curvilinear coordinates. The working fluid is assumed to be a perfect gas and the usual constitutive relations 


3 of 16 


American Institute of Aeronautics and Astronautics 


for a Newtonian fluid are used: the viscous stress tensor is linearly related to the rate-of-strain tensor, and 
the heat flux vector is linearly related to the temperature gradient through the Fourier’s law. A 7th-order 
weighted essentially non-oscillatory (WENO) scheme 26,27 is used to compute the convective flux terms. For 
the viscous flux terms, a 4th order central difference scheme is used. The 3rd order low storage Runge-Kutta 
scheme by Williamson 28 is used for time integration. 

The turbulent inflow is generated using the recycling/rescaling method developed by Xu and Martin 29 
with the recycling station set at 56.7 5% downstream of the inlet to allow for adequate decorrelation be- 
tween the recycling station and the inflow. Dynamical translation operations 30 are applied to the recycled 
turbulence plane at randomly-distributed time intervals to improve the low-frequency characteristics of re- 
cycling/rescaling inflow turbulence generation. On the wall, no-slip conditions are applied for the three 
velocity components and an isothermal condition is used for the temperature with T w ss 0.8 T r . At the 
top and outlet boundaries, unsteady non-reflecting boundary conditions based on Thompson 31 are imposed. 
Periodic boundary conditions are used in the spanwise direction. 

For the current spatial simulation, the boundary layer grows slowly in the streamwise direction, with 
corresponding variations in S from 13.8 mm at the inlet to 24.5 mm near the exit and Reg from 5400 to 
9750. The stremwise computational domain is large enough for the memory of the inflow generation to fade 
out and a nearly uniform acoustic radiation field to be established. In particular, Figure 2 shows that the 
pressure fluctuations at the wall and in the freestream have become nearly homogeneous in the streamwise 
direction after x/Si « 35. 

In the following section, averages are first calculated over a streamwise window of [x a — 0.9A* , x a + 0.9<5i], 
x a = 55<5i, and spanwise locations for each instantaneous flow field; then, an ensemble average over flow 
field snapshots spanning a time interval of approximately 240(5, /[Too is calculated. To monitor the statistical 
convergence, flow statistics are computed by averaging over the whole or half the number of the flow-field 
snapshots and negligible difference (< 1%) is observed between the two. 

Power spectra are calculated using the Welch method 32 with 8 segments and 50% overlap. A Hamming 
window is used for weighting the data prior to the fast Fourier transform (FFT) processing. The sampling 
frequency is approximately 63Uoo/Si, or 4 MHz, and the length of an individual segment is approximately 
53.2(5i/C/oo. Negligible difference in spectral estimation is observed within the reported frequency range, 
when the overall time record is subdivided in 12 segments instead of the baseline number of 8 segments. In 
addition, the DNS methodology has been extensively validated in previous work for supersonic/hypersonic 
turbulent boundary layers 22-25 and for supersonic shock wave-turbulent boundary layer interactions. 33-35 

III. Statistical and spectral analysis 

The main parameters of the boundary layer at the station selected for the statistical analysis ( x a = 55<5^ ) 
are listed in Table 1. The van-Driest transformed mean velocity profile from the DNS is shown in Figure 3. 
The mean velocity conforms well to the incompressible law-of-the-wall and shows a logarithmic region upon 
van Driest transformation. 

A. Characteristics of freestream fluctuations 

In this section, the characteristics of freestream fluctuations are analyzed using the theory of modal analysis, 
which was initially proposed by Kovasznay 36 and further developed by Logan 37 and Smits & Dussauge. 38 
According to Kovasznay, 36 the fluctuations at any point are represented as a superposition of three modes: 
the sound-wave mode (also referred to as the acoustic mode), the entropy mode (also called the entropy 
spottiness or temperature spottiness mode), and the vorticity mode. Modal analysis has been successfully 
applied by Donaldson & Coulter 9 to evaluate the level and the sources of disturbances in AEDC Tunnels A 
and B, and by Masutti et al. 13 to characterize the freestream fluctuations in the Mach 6 Hypersonic Wind 
Tunnel H3 of the von Karman Institute. 

Tables 2 and 3 list the freestream values of several fluctuating flow variables for the Mach 6 and 2.5 DNS, 
respectively. The normalized velocity fluctuations at Mach 6 are larger than those at Mach 2.5. Yet, the 
rms fluctuations in either velocity component are less than approximately 0.2%. Therefore, it is conceivable 
that the dynamics of acoustic fluctuations in the freestream is still linear in nature. However, further work 
is necessary to establish the linearity of the acoustic field. The fluctuations in thermodynamic variables 
are stronger than the velocity fluctuations and also increases from Mach 2.5 to Mach 6. At Mach 6, the 
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rms pressure fluctuations are approximately 2% of the mean pressure value. Moreover, it is shown that 
S rms/R << Prms/P > an d by using the values listed in Table 2, the following isentropic relations are satisfied 


T 

Prms 
P 

indicating that the freestream fluctuations are nearly isentropic and the contribution from the entropy mode is 
minimal. The relative importance of acoustic mode and vorticity mode in the freestream can be demonstrated 
by the ratio of dilatation variance ( dui/dxi )' 2 and vorticity variance D'S2', which are representative of the 
acoustic and vorticity mode, respectively. The large value of (dui/dxi )' 2 implies the dominance of 
acoustic mode over the vorticity mode in the freestream. 

The dominance of the acoustic mode over the other two modes confirms that a purely acoustic field in the 
freestream is successfully isolated by the present DNS. In typical high-speed (noisy) wind tunnels, however, 
multiple disturbance sources exist and all three modes contribute to the freestream fluctuations. 9, 13 The 
present simulation thus provides a unique opportunity to study the generic spectral features and production 
mechanisms particular to the acoustic disturbance. 


7 ~ 1 Pr, 
7 % 

1 z/ 

Jr rms 

7 P 


(1) 

( 2 ) 


B. Intensity of pressure fluctuations 

Figure 4a shows that the wall pressure rms normalized by the local wall shear for the Mach 6 case is about 
2.8, which is close to the corresponding value of approximately 2.6 at Mach 2. 5. 21 The normalized value is 
also close to the value p' Wirms /T w = 3 based on the model by Bies 39 and those given by the DNS of Guarini et 
al., 40 Spalart, 41 and Bernardini & Pirozzoli 42 at lower Mach numbers. Outside the boundary layer, p' rms /T w 
approaches a constant value of about 0.9 for the Mach 6 case, which is significantly larger than the Mach 
2.5 case with p' rms ^ /t w ~ 0.4. The variation of freestream value with Mach number is consistent with the 
trend predicted by the experimental data reported by Laufer 1 (Fig. 5). 

Figure 4b shows the rms pressure profile normalized by the freestream dynamic pressure (<?oo = 1/2 PooU^ = 
1/27 PooM^q). While Prms/qoo shows a strong dependence on Mach number within the boundary layer, it 
nearly collapses the data in the freestream value with p(, mSi00 /<7oo ~ 0.25%, suggesting the potential scaling 
of p' rms oo/poo oc M over the Mach number range covered by these two simulations. However, additional 
data points are necessary to establish the validity of such simple scaling. 

C. Frequency spectra 

The frequency spectrum of the pressure fluctuations is defined as 

1 r°° . 

= — y p'(x,y,z,t)p'{x,y,z,t + r)e luJT dr (3) 

Figures 6a and 6b show the pressure spectrum at the wall and in the freestream, respectively. The 
pressure spectrum is normalized so that the area under each curve is equal to unity. For reference, straight 
lines with slopes of 2, —1, —7/3, and —6 are also included to gauge the rate of spectral roll-off across 
relatively low, mid, and high frequencies, respectively. The wall-pressure spectrum for both Mach numbers 
varies rather weakly with frequency as u> — ► 0. The absence of the more rapid and incompressible u> 2 scaling 
at low frequencies in the wall-pressure spectrum is consistent with the measurements by Beresh et al. 19 and 
the DNS by Bernardini & Pirozzoli 42 at supersonic Mach numbers. At high frequencies, the wall spectrum 
for the Mach 6 case has significantly higher energy than the Mach 2.5 case, and the spectrum for both 
cases exhibits a slightly more rapid decay than the u>~ 5 scaling predicted theoretically by Blake. 43 Unlike 
the Mach 2.5 case, the wall-pressure spectrum for the Mach 6 case shows an apparent u>~ n behavior (with 
n ss 0.7 — 1) corresponding to a universal overlap region. 44 Similar behavior is not obvious for the Mach 2.5 
case and the DNS results (M^ = 2, 3, and 4) reported by Bernardini & Pirozzoli. 42 

Similar to the wall pressure, the freestream pressure fluctuations for both Mach number cases weakly 
decay as ui — ► 0 and exhibit an approximately w -5 roll-off at high frequencies. The freestream spectrum for 
the Mach 6 case has significantly higher energy than that for the lower Mach number case at u>5/U 0 a > 3. 
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While the freestream spectrum for the Mach 2.5 DNS and the experiments by Masutti et al. 13 has an 
observable region of of slope close to —7/3, a similar region is less evident for the present Mach 6 DNS. 

To illustrate the distribution of energy among various frequencies, Fig. 7a shows the pre-multiplied 
pressure spectra at selected heights above the surface. It is shown that the pressure spectra in the inner layer 
(including the wall, buffer layer, and log layer) have a dominant hump centered on luS /U x « 8 (or fS/Uoo 
1), which is the characteristic frequency of the energetic vortical structures within the boundary layer. As 
one moves away from the wall into the outer layer, the peak gradually shifts to lower frequencies as spatial 
intermittency becomes more important. In the freestream, where the the pressure signal is predominantly 
acoustic, the peak of the spectrum is centered at a frequency of ujS/Uoo ~ 3 (i.e. / « 108 KHz), indicating 
that the characteristic frequency of the acoustic mode is significantly lower than that of the vorticity mode. 
Similar variation in pre-multiplied pressure spectrum with wall-normal distance is observed for the Mach 2.5 
case. Figure 7b further compares the pre-multiplied spectra for the two Mach number cases at the wall and 
in the freestream. While the wall spectrum is centered on nearly the same frequency ujS/Uoc 8 at both 
Mach numbers, the freestream spectrum for the Mach 6 case peaks at significantly higher frequency than the 
Mach 2.5 case. The reduced gap in the dominant frequency between pressure signals at the wall and in the 
freestream can be explained by the ‘eddy Mach wave radiation’ concept 45,46 as discussed in Subsection D. 


D. Space-time correlations 


The statistical properties of the pressure field are investigated through the space-time correlation coefficient 
defined as 


C pp ( Ax, Ay , At) 


p'(x, y , z, t)p'( x + Ax, y + Ay, z,t + At) 
(p ,2 {x,y,z,t)J 7 (p' 2 {x + Ax,y + Ay,z,t + At)) 


1/2 


(4) 


where Ax and Ay are spatial separations in the streamwise and spanwise directions, respectively, and At is 
the time delay. 

Figures 8a and 8b display the contours of streamwise-spanwise correlation C pp (Ax, Ay, 0) of the pres- 
sure fluctuations at the wall and in the freestream, respectively, for the Mach 6 case. The contours of 
C pp (Ax,Ay,0) for the Mach 2.5 DNS 21 (dashed lines) are also included to highlight the Mach number de- 
pendence. At both heights, the contours are approximately circular for small spatial separation but become 
elongated in the spanwise direction for large separation distances, indicating that the small-scale pressure- 
carrying eddies are nearly isotropic while the large-scale eddies become more coherent in the spanwise 
direction. Relative to the wall pressure fluctuations, which are primarily vortical, the acoustic radiation in 
the freestream has a significantly greater extent in both in-plane directions, indicating a shift in characteristic 
length scales between those pressure signals. At the wall, the contours exhibit minor Mach number depen- 
dence when nondimensionalized by the boundary layer thickness. In the freestream, however, the large-scale 
pressure-carrying eddies for the Mach 6 case become less coherent and more isotropic compared to the lower 
Mach number case. 

The space-time correlation contours C pp (Ax, 0, At) of the surface and freestream pressure fluctuations 
are shown in Fig. 9a and Fig. 9b, respectively. The shape of the contours at both locations indicates 
the convective nature of the pressure field, which is characterized by downstream propagation of coherent 
pressure-carrying eddies. The overall larger inclination of the space-time correlation contours for the Mach 6 
case indicates the pressure-carrying eddies propagate with a larger convective velocity relative to the Mach 
2.5 case, at least at low to intermediate frequencies. 

From the standpoint of transition analysis, the convection velocity is an important characteristic of the 
stochastic acoustic field in the freestream, since the receptivity characteristics are known to be sensitive to 
the orientation of the plane wave disturbance. 15 To further demonstrate the convective nature of the pressure 
signal as well as its scale dependence, the convection speed is computed using the space-time correlation 
coefficient C pp { Ax, 0, At) following a similar procedure as Bernadini and Pirozzoli. 42 For a given time delay 
At (corresponding to a frequency u> = 2n/ At), the convection speed U c is defined as the ratio Ax\/ At taken 
at the value of Aaq where a local maximum of space-time correlation coefficient is attained. Figure 10a 
plots the convection speed as a function of frequency at several selected heights for the Mach 6 DNS. The 
convection velocity shows a frequency dependence, especially close to the wall. At the wall (z = 0), in the 
viscous sublayer (z + = 5), and in the buffer layer (z + = 17), the convection speed is U c ~ 0.8f7oo for the 
large-scale disturbances (associated with lower frequencies) and U c ~ OTf/oo for the small-scale disturbances 
(associated with higher frequencies), which is significantly larger than the local mean velocity. The values 
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of convection speed at the wall are consistent with those widely quoted for low-speed flows 47-49 and those 
reported by Bernadini and Pirozzoli 42 for turbulent boundary layers at Mach 2, 3, and 4. In the log layer 
(z/8 = 0.16) and the outer layer (z/8 = 0.75) of the boundary layer, the convection velocity becomes less 
frequency-dependent and is nearly equal to the local mean velocity value. In the freestream (z/S = 3.2), the 
pressure-carrying eddies propagate at significantly smaller speeds with U c s=s 0. 63Uoo. A similar variation of 
convection speed with the wall-normal distance is observed in the Mach 2.5 case. 21 In addition, Fig. 10b 
further compares the convection speeds at the wall and in the freestream between the Mach 6 and Mach 2.5 
cases. It is seen that the convection speed both at the wall and in the freestream is larger for the Mach 6 
case, consistent with Fig. 9a and Fig. 9b. 

To represent the overall convection speed of the pressure-carrying eddies, we define the bulk convection 
speed (14) as that associated with the frequency where the pre-multiplied pressure spectra (Fig. 7a) attain 
their maximum. The bulk convection speed for the Mach 6 case is Ub ~ 0.71L4o at the wall and 14 ~ 0.63C4o 
in the freestream, compared with Ub ~ 0.66 at the wall and Ub ~ 0.41/oc in the freestream for the Mach 
2.5 case. Figures 11a and lib show the variation of bulk convection speed with freestream Mach number at 
the wall and in the freestream, respectively. While the data by Kistler & Chen 17 shows that the convection 
speed at the wall plateaus for > 2, the current DNS and the data reported by Bernadini and Pirozzoli 42 
show a weak increase with the freestream Mach number. Both the current DNS and the experimental 
measurements by Laufer 1 show that the bulk convection speed in the freestream increases with freestream 
Mach number, while the value reported by Kendall 50 at = 4.5 is relatively higher. As expected, all the 
freestream bulk convection speeds fall within the region where M r > 1, with M r = (L4o — Ub)/cioo- 

Figures 12a and 12b show that the instantaneous pressure field in the freestream exhibits a ‘Mach-wave’ 
radiation pattern with a preferred direction. The fact that the radiation wave front for the Mach 6 case is 
shallower than the Mach 2.5 case is consistent with the experimentally measured trend. 1 By assuming that 
the acoustic radiation consists of plane waves generated by ‘effective’ acoustic sources within the boundary 
layer, the convection speed of acoustic sources, U 8 , can be calculated. It is shown that the angle between 
the flow direction and plane-wave normal direction, 9 n can be determined by 1 


cosd n 


—jMoo 


u'/Uqq 

P'/Poo 


-"fMoo 


V 

rms 

Prms 


/Uoo 

/ Pea 


(5) 


where the negative sign specifies that the plane wave is radiating away from the boundary layer with 9 n > 90° . 
Once 9 n or the wave-front orientation angle 9 (shown in Fig. 12) is known, the velocity U s of sound sources 
that produce the waves can be determined by the ‘Mach angle’ relation 


1 

sin9 


= M r 


Uoo -Us 

O'oo 


( 6 ) 


By using the values of u' rrrls /U x and p' rrns /Poo listed in Table 2 for the Mach 6 case, Eq.(5) gives a value 
of 9 n « 120° (correspondingly 9 = 30°), which agrees with the approximate inclination of the wave fronts 
from the numerical Schlieren image in Fig. 12a. Moreover, Eq. (6) gives the relative Mach number M r = 2 
with the corresponding convection velocity of ‘effective’ radiation sources U s ~ 0.66L4o, which is slightly 
greater than the bulk convection speed Ub = 0.63 calculated based on the space-time correlation (Fig. 11). 
An approximate agreement between the two convection speeds was also found for the Mach 2.5 case 21 with 
the ‘effective’ source speed U s « 0.4 and M r = 1.5. 

The supersonic value of relative Mach number for the ‘effective’ radiation sources as well as the match 
between the source convection velocity with the bulk convection speed Ub supports the basic concept of ‘eddy 
Mach wave’ 45,46 and shows that the ‘Mach wave type’ radiation is produced by eddies which are convected 
supersonically with respect to the freestream. Although a quantitative validation of the theory is beyond 
the scope of this paper, the theory can be used to explain the observed variation in radiation intensity, 
directionality, and convection speed with freestreem Mach number. At low supersonic freestream Mach 
numbers, sources that contribute to the radiation field are primarily slowly moving ones, the convection 
velocities of which are supersonic relative to the freestream. As the Mach number increases, additional 
faster moving turbulent eddies acquire supersonic relative speeds and start to take part in the radiation 
process, acounting for the increased ‘effective’ source convection speed and enhanced radiation intensity. 
This explains the increase in p' rmSt00 /T w from M w = 2.5 to = 6 as found in Fig. 4a. Moreover, the 
increased fraction of the inner layer that can radiate to the free stream contributes to a reduced gap between 
the peak frequency of fluctuations near the surface and within the free stream as shown by Fig. 7b. 
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IV. Conclusions 


DNS is used to examine the pressure fluctuations generated by a hypersonic turbulent boundary layer 
with a nominal freestream Mach number of 6. Detailed analyses of the statistical characteristics of pressure 
fluctuations, including the fluctuation intensities, frequency spectra, space-time correlations, and convection 
velocities, have been carried out. The DNS results are compared with the recently reported Mach 2.5 DNS 21 
to highlight the variation of acoustic characteristics with Mach number. It is found that the Mach 6 DNS 
exhibits increased radiation intensity, enhanced energy content at high frequencies, shallower orientation 
of wave fronts with respect to the flow direction, and larger convection velocity relative to the Mach 2.5 
case. These variations in the freestream pressure field with Mach number agree well with the experimentally 
measured trends 1 and are consistent with the ‘Mach wave radiation’ concept. 45,46 
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Table 1. Boundary layer properties at the station selected for the analysis ( x a = 55<5i) of the 
acoustic field for the current Mach 6 DNS and Mach 2.5 DNS. 21 


Moo 

E/oo(m/s) 

Poo (kg/ m3 ) 

Too(K) 

T w ( K) 

R eg 

Re T 

r £82 

0(mm) 

H 

6 (mm) 

Si (mm) 

5.85 

870.2 

0.039 

55.0 

300.1 

9659.2 

464.2 

1783.3 

0.968 

13.6 

24.4 

13.8 

2.5 

823.6 

0.1 

270.0 

568.0 

2834.8 

509.9 

1656.9 

0.583 

4.14 

7.69 

4.0 


Table 2. The disturbance field in the freestream for the Mach 6 DNS at z/5 = 3.2, Re T — 464. 


^ rms 17 rms Prmsl P Prmsl P r ^'rms/ r ^- 

1.26 x 1CT 3 0.99 x 10" 3 2.00 x 10" 3 2.06 x 10" 2 1.47 x 10" 2 5.86 x 10" 3 

(Pu)rms/W T lrms/ T t P't,rms/Pt S 'rms/ R 

1.38 X 10" 2 1.76 x 10 -3 6.15 x 1CT 3 1131.2 1.75 x 10~ 4 


Table 3. The disturbance field in the freestream for the Mach 2.5 DNS at z/5 = 2.8, Re T = 510. 21 


^ rms ^ rms ^ rms / ^ Prms/P Prmsl P ^rms 1^ 

7.56 x 10 -4 4.90 x 10" 4 8.09 x 10" 4 3.96 x 10" 3 2.83 x 10" 3 1.13 x 10" 3 

(P u )'rms/W Tl rrns /T t Pt,rms/Pt ^o ^' 2 S 'rms/ R 
2.32 x 1CT 3 5.75 x 10“ 4 2.01 x 10" 3 1797.8 2.19 x 10" 4 
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Figure 1. Computational domain and simulation setup. The reference length 5i is the thickness of 
the boundary layer (based on 99% of the freestream velocity) at the inflow plane. An instantaneous 
flow field is shown in the domain, visualized by an iso-surface of the magnitude of density gradient, 
\Vp\Si/ P ao = 0.98, colored by the streamwise velocity component (with levels from 0 to (Too, blue 
to red). Note that x, y, and z are, respectively, the streamwise, spanwise, and wall-normal 
coordinates. 



Figure 2. The streamwise distribution of the normalized pressure fluctuation rms at the wall 
(z/5 = 0) and in the freestream ( z/5 = 3.2) for the Mach 6 DNS. 
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Figure 3. van Driest transformed mean velocity profile for the Mach 6 DNS (k = 0.41, C = 5.5). 




Figure 4. Pressure fluctuation rms profile for the Mach 6 and Mach 2.5 DNS. 21 (a) normalized 
by wall shear stress r w j (b) normalized by freestream dynamic pressure q a Q . 
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Figure 5. Intensity of freestream pressure fluctuation for the Mach 6 and Mach 2.5 DNS 21 com- 
pared with the experiments by Laufer. 1 




(a) Wall (b) Freestream 

Figure 6. Normalized frequency spectrum of computed pressure signal at selected heights for the 
Mach 6 and Mach 2.5 DNS. 21 The freestream are taken at z/5 = 3.2 for the Mach 6 DNS and 
z/8 = 2.8 for the Mach 2.5. 
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Figure 7. Pre-multiplied power spectrum of pressure signals, (a) at selected heights for the Mach 
6 DNS; (b) comparison between Mach 6 and Mach 2. 5. 21 The pressure spectrum is normalized 
so that the area under each curve is equal to unity. The freestream are taken at z/5 = 3.2 for the 
Mach 6 DNS and z/5 = 2.8 for the Mach 2.5. 



Ax/§ Ax/§ 


(a) Wall (b) Freestream 

Figure 8. Streamwise-spanwise correlation coefficient of the pressure signal at selected heights 
for the Mach 6 DNS (Solid line) and Mach 2.5 DNS 21 (Dashed line). Contour levels vary from 
0.1 to 0.9 with increments of 0.1. The freestream are taken at z/5 = 3.2 for the Mach 6 DNS and 
z/5 = 2.8 for the Mach 2.5. 
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Figure 9. Space-time correlation coefficient of the pressure signal at selected heights for the Mach 
6 DNS (Solid line) and Mach 2.5 DNS 21 (Dashed line). Contour levels vary from 0.1 to 0.9 with 
increments of 0.1. The freestream are taken at z/5 = 3.2 for the Mach 6 DNS and z/5 = 2.8 for the 
Mach 2.5. 




Figure 10. Convection speed of the pressure fluctuation as a function of frequency at selected 
heights, (a) at selected heights for the Mach 6 DNS; (b) comparison between Mach 6 and Mach 
2. 5. 21 The freestream are taken at z/5 = 3.2 for the Mach 6 DNS and z/5 = 2.8 for the Mach 2.5. 
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(a) (b) 

Figure 11. Bulk convection speeds of the pressure fluctuation as a function of freestream Mach 
number, (a) at the wall; (b) in the freestream. Symbols: squares, Kistler & Chen; 17 deltas, 
Bernardini &: Pirozzoli; 42 right triangles, Laufer; 1 left triangles, Kendall; 50 diamonds, the present 
DNS and Mach 2.5 DNS. 21 Lines: M r = 1. 



Figure 12. Numerical Schlieren image based on instantaneous flow field, (a) Mach 6 (b) Mach 
2. 5. 21 Contour levels are selected to emphasize disturbances in the freestream. In (a), the dashdot 
line indicates the orientation of the radiation wave front; n is the direction normal to the wave 
front; U is the flow direction; 9 n is the angle between U and n. 
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